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Explicit computation of Drinfeld associator in the case of 
the fundamental representation of gl(N) 



P.Dunin-Barkowski* A.Sleptsovj A.Smirnov^ 



ABSTRACT 



We solve the regularized Knizhnik-Zamolodchikov equation and find an explicit expression for the Drinfeld associator. We restrict to 
the case of the fundamental representation oi gl{N). Several tests of the results are presented. It can be explicitly seen that components 
of this solution for the associator coincide with certain components of WZW conformal block for primary fields. We introduce the 
symmetrized version of the Drinfeld associator by dropping the odd terms. The symmetrized associator gives the same knot invariants, 
but has a simpler structure and is fully characterized by one symmetric function which we call the Drinfeld prepotential. 



1 Introduction 

It is well known that the polynomial invariants of knots can be evaluated with the help of Chern-Simons theory 
[I]. We recall some essential facts and constructions. 

The Chern-Simons action is given by Scs = / t^R {A A dA + | A A A A A). Chern-Simons theory is a three- 
dimensional topological quantum field theory with gauge group G; traces are taken in a fixed representation R 
of this group. Throughout this paper we let G = GL{N). One can consider the gauge invariant Wilson loop 
operators, which are defined as follows: 



The vacuum expectation value ([T]) does not depend on the realization C of the knot in M^, but only on the 
equivalence class of the knot K therefore, ( W{K) ) defines a knot invariant. Then in order to evaluate the integral 
([I]) one fixes some particular gauge. There are two relevant gauges: the holomorphic gauge [213] and the temporal 
gauge I4j 5^ . 

If the temporal gauge is chosen then integral ([T]) gives a polynomial invariant, namely it is proportional to the 
HOMFLY polynomial which specializes to the Jones polynomial if one fixes a dimension N — 2 gauge group. In the 
holomorphic gauge case, the vacuum expectation value ^ is equal to the celebrated Kontsevich integral. Due to 
the gauge invariance of Chern-Simons theory these answers should coincide, i.e. the Kontsevich integral computed 
for a particular knot (or link) is equal to the corresponding HOMFLY polynomial of this knot (or link). However 
the Kontsevich integral is given by a complicated infinite power series ([6l|3]). It is impossible to sum this series in 
general and thus it is difficult to check this explicitly. 

On the other hand the Kontsevich integral for knots and links can be calculated as a trace (over a representation 
R of the gauge group G) of a certain product of two basic operators: an R-matrix and a Drinfeld associator <I> 
[3 E]- The R-matrix has a very simple structure, while the associator looks almost as complicated as the 
Kontsevich integral itself. It is given by an infinite power series in two non-commuting variables A and B with 
multiple zeta values as coefficients. In order to compute the Kontsevich integral in this way one has to sum the 
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series for the associator. This paper addresses the problem of summing the associator series for the simplest case 
of the fundamental representation (see also |H1E])- The key is the Knizhnik-Zamolodchikov equation, which is an 
additional constraint satisfied by the correlation function of the conformal field theory associated with an affine Lie 
algebra at a fixed level (see ^0|). The solutions to this equation are the correlation functions of primary fields (see 
[T2]V However, in practice they are rather difficult to solve directly, except for the four-point functions, which are 
directly related to associators [B]. The only difference between associators and these four-point functions is that 
the associator is a solution of a regularized version of the equation. In addition, there are actually three slightly 
different versions of associators and their corresponding equations. 

The problem of explicitly evaluating the Kontsevich integral may seem to be just a technical problem. However 
we emphasize that the search of explicit expressions for the associator $ is very important and is a conceptual 
problem. The point is that the associator plays a key role in many areas of mathematics and theoretical physics 
such as Hopf algebras, conformal field theory, quantum groups, braid theory and many others. 

In this paper we present an explicit expression for all forms of the associator in the case of the fundamental 
representation of gl{N). The coefficients of one of these forms turn out to coincide with certain components of the 
WZW conformal block for primary fields. We apply the obtained results to a few simple knots to check that the 
expression for the Wilson loop in the two gauges mentioned above actually coincide. 

The odd part of the associator is represented by the values of multiple zeta functions at some odd points, 
and therefore does not enter the expression for knot invariants. We show that dropping the odd terms leads to a 
symmetrized version of the Drinfeld associator whose structure is simpler: it is fully characterized by one function 
even in Planck constant. We call it the Drinfeld prepotential. 

Our paper is organized as follows. In Section [2] we introduce the Knizhnik-Zamolodchikov equation and its 
regularized form, and then we list three different forms of this equation needed for the knot theory. In Section [3] 
we provide explicit solutions to all forms of the equation from the previous section, for the case of the fundamental 
representation of gl{N). In Section |4] we apply the obtained results to several particular examples of knots and 
find that they lead to invariants coinciding with the HOMFLY polynomials. In the final section, we introduce the 
symmetrized version of the associator which does not contain odd zeta function values. 



2 Regularized Knizhnik-Zamolodchikov equation for associator 

Associators appearing in combinatorial construction of Kontsevich integral coefficients are related to the so-called 
Knizhnik-Zamolodchikov equation. This is an equation coming from conformal field theory \n . Consider the 
four-point function in WZW model: 

G{z) = {g{0)g{l)g{z)g{(xi))w zw 
Knizhnik-Zamolodchikov equation is the Ward identity for this four-point function: 

dG{z) ^ B 



dz 

where A and B have the form (we assume summation over repetitive indices): 

yi = r"«)(r")*«)i, B ^ 1 ® T'' ® {T'')* 

G{z) is a rank-six tensor of the same type as A and B. Also note that everywhere where we write a product of 
two such rank-six tensors, we assume the following: 

{ABTs^;: ■■= a-p;b^,:^ (3) 

The perturbative approach for solving of the equation ^ is discussed in Appendix |Aj The solution in the first 
two orders of perturbative expansion is also given. Solutions of slj^ Knizhnik-Zamolodchikov equations at level 
zero are given in |14j . However we are interested in general situation and in regularized Knizhnik-Zamolodchikov 
equation. Moreover we need explicit expression for assocaitor to apply the result to the knot theory 

This solution G{z) diverges at points and 1, so we make the following transformation to eliminate the 
divergences: 

G{z)^{l~ z)-^^ <^{z)z^'^. (4) 
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The regularized equation has then the form 



dz z z 



(5) 



The solution of this equation taken at point z = 1 gives us the associator. In Section |3] equation ([s]) is solved in 
the case of fundamental representation of gl{N). 



2.1 T ypes of associators 

Actually, to construct the coefficients of the Kontsevich integral combinatorially (see 013]), one needs the i?-matrix 
and three types of associators. These three types of associators correspond to three types of orientation of strands 
of a braid and lead to three slightly different forms of equation They are listed in the following table: 



Orientation 


KZ- 


- equation 




Transformation 


Associator equation 


a) 


ill 


dGiz)^, 
dz V 


A B \ 
z 1 — z J 


G{z) 


G(z) = (1 - zf^ $(z)z'"'^ 


dz z z 


in 


dG{z) ^ 
dz 


^A B 

\ z 1 — z 


)g(z) 


G(z) = (1 - zy"^ 


^ = Z^il-^rA{l-zy'^<^>+^^A 
dz z z 




III 


dGiz)^, 
dz \ 


A B \ 

z 1 — z/ 


G(2) 


G(z) = (1 - zy"^ $(0)z'^^ 


— = - (1 - zf" ^ (1 - ^) $ - - $A 
dz z z 



(6) 



3 Solutions 

In this section we solve the corresponding equation for each of three cases ([6]). We restrict our considerations to 
the case of fundamental representation of gl{N) algebra only. 

3.1 Type a 

Let us consider type a from table [6j In this case A and B take the following form 

A^T"® (T")*® 1, B = 1 «)T° (g) (T°)* , (7) 

where, since we are considering the fundamental representation of gl{N), 

T^=g,,, [T^y =g^„ (8) 

where index a corresponds to a pair of indices i and j each of which goes from 1 to and gij is a matrix with 
components equal to except the (i,j)-th one, which is equal to 1. Then we get the following relations for A 
and B: 

A^ = B^ = {ABf = {BAf = 1. (9) 
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In the paper [5] it is proven that associator $ is given by the infinite series in A and B ( 48 1 . However there are 
relations ([9]) on ^ and -B, therefore associator can be written as finite polynomial: 



$ = $1 + $2^ + + $4^5 + ^^BA + ^qABA, 



(10) 



where $i are some coefficients and do not depend on A and B. Now let us substitute (10) in equation ([5|, then it 
turns into a system of differential equations, which can be represented in the matrix form: 



where 



Matrix M has the form: 



where 



dz 



M(z)$, 



$2 
$3 
$4 
$5 



(11) 



(12) 





( 




s2 





— cs 


cs 


— s^ 













cs 


-S2 





— cs 






M = 





—cs 







-C2 


cs 




(13) 




cs 









— cs 


-1 








—cs 







CS 





C2 








V ' 


cs 


— CS 


-1 


c^ 


° J 






c := cosh {h\og 


(1-^)) 


1 

" 2 


((- 


z)V 


(1-z)- 


-), 


(14) 


s :— sinh {hlog 


(1-^)) 


1 

~ 2 


((- 




(1-z)- 


-')■ 


(15) 



In order to solve equation (11 1 one can change the basis in the following way: 







: 5$, 








1 








1 


1 








1 


1 








1 


-73 




73 





73 





-1 


-1 


-1 


2 


-1 


2 


V3 


-V3 


73 





-73 





-1 


1 


1 


2 


-1 


-2 



(16) 



(17) 



Then matrix M takes the form 

/ 



M = SMS-^ 









^[l-z] 



-2h 











{l-zf 

3 
2 






f (1-z) 



2h 



3 
2 



Hence, matrix equation (11 1 turns into the foUowing system of differential equations: 



dz 



|;i$4 + #/i(i-^)""'$3, 



with the following initial conditions ($(0) = 1): 



5^(0) = -1, 
5^i(o) = -1 



From (19) and (20l one can obtain conditions for derivatives of $i(0) using $(0) = hM{Q)<^{Q) + fiM(0)$(0) 



$3(0) 



1 + 2/1' 

3/12 



Thus systems of equations ( 19]|2T ) can be reduced to a system of hypergeometric equations: 

z (1 - Jij + (1 + 2fi - z) 5^3 + /i^Jjj = 0, 

z {1 - z)¥^ + {1 + 2h - {1 + ih) z)¥i - = 
z (1 - z) 5^ + (1 - 2/1 - z) $^ + h^¥5 = 0, 



(22) 



z (1 - z) $6 + (1 - 2/i - (1 - iH) z) $6 - 3;i2$6 = 0, 
We get a system of differntial equations of the second order. Solutions of this system considered with initial 



conditions (20 21 ) are equivalent to solutions of ( 19 ). Solutions are given by the following hypergeometric functions: 



$1 = 1, 
5^ = 0, 

^3 = -VI 2Fi{K-h,l + 2h:z) , 
^A = - 2F1 {h, 1 + 2h; z) , 
5J = \/3 2i^i (fi, -n, 1 - 2S; z) , 
^6 = - 2F1 i-h, -3h, 1 - 2h; z) 



(23) 



We have checked with the help of a computer algebra system that solutions (23) indeed obey original equations 



(19) with initial conditions (20) 



Thus we see that coefficients of associator of this type are given by linear combinations of values of hypergeo- 
metric functions taken at the point z = 1, which can be rewritten in terms of Gamma function and trigonometric 
functions. Note that the answer in this case does not depend on order N of algebra gl{N): 



$1(1) = 1, 

5^(1) = 0, 

¥3(1) ^ -V3 2F1 {h, -H, 1 + 2h; 1) 
5^(1) = - 2Fi{h,3h,l + 2h;l), 
¥5il)^VS2Fi{h,-h,l-2h;l), 

M^) = - 2Fii-h,-3n,i-2n;i) 



$1(1) = 1, 



$2(1) 
53(1) 
$^(1) 

5^(1) 
5^(1) 



= 0, 



V3 



r(i + 2fi)2 



r(i + fi)r(i + 3fi)' 



(24) 



cos(7r/i) ' 



r(i - 2hf 



r(i-fi)r(i-3fi)' 
1 



cos(7r/i) 
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Thus we obtained explicit expressions for $. To get expression for associator one needs to take ^ — S and 



then use ( 10 1 



3.2 Type b 

Now we briefly present results for type b from table |6] In this case we have: 

A = Qij (gjgij (E)l, B = 1(E) g^j (g) 
Relations for A and B are the following ones 

A^ = NA, 
B^ = NB, 
BAB = B, 
ABA = A 
Then the associator takes the form 

$ = $1 + + + •^iAB + $5^^. 
By analogy with the previous case we write 



$ 5$, 



for 



S- 



1 














1 


N 





1 





1 








-1 





1 





N 


1 





1 

Ar2 


1 

N 


1 

N 


1 


1 



Then we have 
5^ = 1, 



$2 = 2F1 {h, -h, 1 - Nh; z) , 

1 



Fi {{N - l)h, {N + l)h, 1 + Nh; z) 



- 1 

51= 2Fi{h,-h,l + Nh;z), 



$5 = ^ 2F1 {-{N - l)h, ~{N + l)h, 1 - Nh; z) 



$1(1) = 1, 

5^(1) = 
^3(1) = 
5^(1) = 



T{l-Nhf 



r(i - Nh- h)T{i - Nh + h)' 

1 Nsm{nh) 



iV2 - 1 sin(7riV;i) ' 

T{l + Nh)^ 

r(i + Nh- h)r{i + Nh + hy 

1 Nsin{-nh) 



7V2 sm{'KNh)' 



(25) 



(26) 



(27) 
(28) 



(29) 



(30) 
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Here we also see that coefRcients of associator of this type are given by hnear combinations of vahics of hyper- 
geometric functions taken at the point z — 1. 

Note that the coefficients of the associator in this case coincide with certain components of the WZW conformal 



block for primary fields in the fundamental representation. Namely, we have 



(-) 



(+) 



1 



$3 up to singular 



factors and substitution h — where F^ ^ and F^^' are the components of the WZW conformal block for primary 

fields taken from Section 15.3.2 of [12]. The other two coefficients $4 and $5 differ from the previous two just by 
substitution N ^ —N. 



Again, to get the expression for associator one needs to take $ = 5* ^<f> and then use (27) 



3.3 Type c 

Associators of type c allow us to calculate the simplest knots like in Figure [l] In this case we have: 



2k+1 I 




Figure 1: (2k+l)-foil knot 



A = cjij (g) gji(g)l, B = l(g) ® g^j 
Relations for A and B are the following ones: 

A^ = 1, 

, B^ ^ NB, 

BAB = B, 

Then the associator takes the form 

$ = $1 + + + ^4AB + ^^BA + ^qABA. 
By analogy with the previous case we write 

$ = S$, 



(31) 



(32) 



(33) 
(34) 
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for 



1 

















-N 


1 














1 


1 


1 


1 


1 


1 


2N 


2N 


2 


2N 


2 


2N 


1 


1 





1 





1 


2{N+1) 


2(Af+l) 


2 


2 


1 


1 


-1 


1 


1 


1 


N 


N 


N 


w 


1 

\^ Af-1 


1 

AT-l 





-1 





1 



(35) 



Then 



$1 = 1, 



$2 = -N, 



$3 



$4 



27V 



2Fiih,-{N -l)h,l + 2h; z) 



1 



$5 = 



2{N + 1 
1 



ai^i (/i, + l)h, 1 + 2/1; z) 



1 



TV- 1 



2F1 {-h, -{N + l)h, 1 - 2h; z) 
2Fi{-h, {N-l)h,l-2h;z) 



$i(l) = 1 

$;(!) = - 

5^(1) 
5^(1) 

5^5(1) 
$6(1) 



iV, 



1 'i''T{h+ l)T{i + Nh) 
2N y/^r{i + Nh + h) ' 

1 4'*r(fi+ i)r(i - A^/i) 
2(7V + i) v^r(i -Nh + h) ' 

1 4-'^r(-fi + i)r(i + 7v;i) 
^iv v^r(i + ivfi-;i) ' 

1 4-'"T(-fi+ i)r(l - A^/i) 



(36) 



/V - 1 V^T{i - Nh- h) 



In this case as in previous ones coefficients of associator of this type are given by hnear combinations of values 
of hypergeomctric functions taken at the point z — 1. 

4 Comparison to known formulas for knots 

In the present section we apply the formulas for the associator obtained in Section [3] to compute Kontsevich 
integral for a few knots. We do this following [5J[T5]. The computations are done in the case of the fundamental 
representation of gl{N). In this case Kontsevich integral is known to coincide with HOMFLY polynomial for 
the given knot, and we use this fact to check the validity of our formulas for the associator. Namely, we check 
that Kontsevich integral computed with the help of these formulas indeed coincides with known expressions for 
HOMFLY polynomial. 

4.1 Hump 

We start with the hump, which is a version of the unknot (for the definition of the hump see [6]). Recall that 
associator is a tensor with 6 indices. Denote the associtaor of type b as $b. Kontsevich integral of the hump 
corresponds, by definition, to the contraction of the indices of $f, corresponding to Figure [2j We will refer to this 
contraction as the trace of the assocaitor. 

Writing this contraction in terms of the indices gives the following expression: 



KI„ 



= tr„ 



Now let us expend the associator as in formula (27) 



(37) 



(38) 
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Figure 2: Hump 



With the help of formula ( 25 ) one can compute traces (in the above mentioned sense) of particular terms of 



(38) in the following way: 



trHump ^ 



trH„„p AB 



1 

N 



trHump = ]^ tr {gi-jgmngnm9ji) = N'^ 



(39) 
(40) 
(41) 



Note that trnu^p here stands for the above mentioned trace of a rank-six tensor, while tr stands for just the ordinary 
trace of a matrix. We also assume summation over repetitive indices. 

Now substituting formulas (30), (29), (39)-(41) to (28) we get the following answer: 

KiHump = 2F1 {~{N - i)n, -(TV + i)n, i-Nn-i)^N 

where we introduce the following notation for a g'-number 



N 



(42) 



— 1 ' 

q-q 



q = e 



This indeed coincides with the known result for the HOMFLY polynomial of the hump (which can be found 
e.g. in [TS]). 

4.2 Series [2,2k+l] 

Consider the so-called torus knots of type [2,2k+l], as shown in Figure [TJ Denote the associator of type c as $c- 
With the help of this associator one can compute Kontsevich integral for knots of this type [51 [TS] : 



where R is in our notation equal to 



0fs 

77' 



(43) 



(44) 



Now expanding the exponential in series, substituting formula (33) for $c and taking into account relations (32), 
one can arrive at the formula 



KI[2.2k+l] = Ki + K2A + K3B + K4AB + + KqABA, 



(45) 



where ni, . . . ,kq are certain expressions in terms of the coefficients $1, . . . , $g from formula (33). 

Then, after taking traces of the individual terms in a way anologous to formulas ( 39|[4T | and substituting 
$1, . . . , <i>6 from formula (32), we get 



KI[2,2k+l] = N 



cos(7r/i(2fc -I- 1)) sin(7r?i) cos{TrhN) + i ■ sm{TrH{2k + 1)) cos(7rS) sin(7r?iA'') 

cos(7rfi.) sm{TrhN) 



(46) 
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Note that this answer is not normahzed with respect to the unknot. Let us introduce q — e^"^^ and normahze with 
respect to the unknot: 

„3 („2k+2 _ -2N+2k , -2N-2k _ -2k-2\ 

if/[2,2.+i]/i^/(0) = '-^ '- (47) 

This expression is in a perfect agreement with known resuhs. For example, it is known that in the AbeUan case 
= 1 the expectation value of the Wilson loop is given by q''^-^\ where 1{K) is the so called linking number of 
the knot K [S]. The linking number can be computed as a sum of the orientations for the crossings of the knot, 
which takes values ±1. In our case - Figure [ij the knot has 2k + 1 crossings with equal orientations, so the linking 



mumber is Z = 2A: + 1. It is an elementary exercise to check that for iV = 1 our formula (47) indeed gives q'^'^'^^. 



In the first nontrivial case N — 2 one can check that the result (47 1 for any choice of fc is a Laurent polynomial 
in q which coincide precisely with the Jones polynomial for the torus knot of type [2, 2fc + 1]. For a general N the 
answer (47) is a polynomial in ^ = q~^ and a Laurent polynomial in q. This two parametric polynomials again 



coincide with known HOMFLY polynomials for these knots as computed for example in [ini [H] . Being a rather 
nontrivial computation, this provides a strong test of our formulas for the associator. 



5 The Drinfeld prepotent ial 

In the previous sections we have computed the Drinfeld associator for the case of the fundamental representation 
explicitly and checked that it reproduces correct HOMFLY polynomials for toric the knots of type [2, 2fc + 1]. In 
fact, as we explain below, the Drinfeld associator contains more information than is used for the knot invariants in 
the following sense: not all components of the Drinfeld associator contribute to the answer for the knot invariants. 

Using the perturbation approach to the solution of the regularized Knizhnik-Zamolodchikov equation (as in[A|, 
T.Le and J. Murakami found the expression for the associator as infinite sum in non-commutative operators A and 
B (without any additional relations of the form i.e. for all representations) with the coefficients in the form of 
the multiple zeta function values j8]: 



Pi k 

*3-l^' + E(2!^) E E (-l)l''lr(pi,gi,...p™,<z™)x 



k=2 



m>0 p>Oq>o 
- |p| + |q| = fc 

i(p) = i(q) = ™ 



(48) 



!(r) = i(B) = „ 
0<r<p, 0<s<c 



Where p = {pi,p2, ■■■,Pm,) is a vector with positive integer components. The length of the vector lijp) — m and 
|p| = SPi- For p and q the with positive integers p > q means that > qi and p > means that > for all 
i. The coefficients t(pi, gi, ...pmj 9m) £^re expressed through the multiple zeta functions as follows: 



r(pi, gi, ...p,„, g„j) = C(l, q\ + 1,1, 1, q2 + 1, + 1) 

Pl-l P2-1 

such that for example r(l, 2) = C(3) and t(2, 1) = C(l, 2). The zeta functions are defined as: 

C(mi,TO2, ...,m„) = ^ fcJ;"lfc2"'"^..fc~™" 

0<fcl<fe2<...<fc„ 

The operators A and B depend on the gauge group and its representation: 



(49) 



(50) 



A^^®\, B = i(»n, 



(51) 



Formula ( 48 ) gives exact answer for the Drinfeld associator in all orders of h, however it is impractical for explicit 



computations of the knot invariants. For such computations one needs to summate exactly series (48) in the 



presence of some relations among A and B given by a representation of the gauge group. In the case considered 
here (fundamental representation of gl{N)) these relations are given by ([£]), (26 ) or (32 ) depending on the associator 
type. The main difficulty is that the values of the multiple zeta functions (50) are difficult to compute (see for 
example [T7]) and there are no chances to summate series (48) exactly in all orders of h by "bare hands" even in 
the presence of additional relations for A and B. 
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Our approach implements the relations from the very beginning: using them, we reduce the KZ equation to 
the system of the hypergeometric equations. The solution of this system at the point z = 1 gives the associator. 



In this way, we were able to summate the series ( 48 ) with nontrivial zeta functional coefficients explicitly and the 



result has the form of simple gamma functional factors ( 24 ) . Using computer we checked explicitly that up to the 



order our result coincide precisely with series coming from ( 48 ) . This should hold for all higher orders as well 



but we do not know how to show it explicitly due to the multiple zeta functions being difficult to work with, as 
mentioned above. 

It is known, that the knot invariants must be rational functions in q — exp{iTrh) for any choice of gauge group 



and its representation. For example our answers for unknot (42) and torus knots of type [2,2fc + 1] (46) are 
obviously of this type. Therefore, the /i-series expansion of knot invariants has the form X]fcLo^'=^'^ where we 
denoted h = itt/i, with rational coefficients Ofc S Q. 

The even values of the multiple zeta functions, i.e. C(™i: 't*2j "m-n) for k = uii +m2, ... + m„-even number, are 
known to have the form ott*^ where a is a rational number a G Q. On the other hand the odd zeta function values 
are always irrational (conjecturally are transcendental Moreover, the odd zeta values are algebraically 

independent over rational numbers, i.e. there is no polynomial relations among this values over Q ( the coefficient 
of polynomials are rational) |17J . 



It is clear from the above formula (48) that the odd (even) zeta values represents the odd (even) coefficients 
of the Drinfeld associator expanded in the Planck's constant h. Therefore, in order to make the coefficients of the 
fi,-series for invariants rational, the odd part of the associator must cancel in the answer. Indeed, the knot invariant 
is given by a trace (over knot projection to two-dimensional plane) of a product of associators and brading matrices 



R as for example in (43). Brading matrix is given by exponent (44) and its coefficients of /i-expansion are rational 
for all representations. Therefore, the contribution of the odd part of associator to the knot invariant is given by 
polynomial combinations of odd zeta values with rational coefficients, and as mentioned above such a combination 
is rational only if all the coefficients of the polynomial are zero, i.e. the coefficients of odd zeta values cancels. 
Calculations up to shows the following: every time we are computing the knot invariant by taking the 



trace of a braiding, for example as in (43), the odd zeta values enter the answer with a group factor (i.e. traces 
words made from A and B) that are exactly zero due to celebrated STU and IHX relations for chord diagrams llj. 
Therefore, for knot theory, we can use only the even part of the associator. We found, that dropping the odd terms 
in some clever way, the expression for associator possess valuable simplification. Indeed, consider the logarithm F 
of the associator 

$ = exp(i^(;i)) (52) 
Introduce then the symmetrized associator (i.e. we drop all odd terms): 

<i>,^e^p[^iFih) + Fi~h))), (53) 

It should be clear from the discussion above , that this "even" associator leads to the same knot invariants, but its 
structure is simpler. For example, in the fundamental representation case, for the associator of type (a) specified 
by the relations ([9|: 

A'^ = B^ = {ABf = {BAf = 1. (54) 

we obtain: 

<i>s=exp(^Ts{h^)[A,B]^ (55) 

with: 

^ ' 2VN^-1 V [1 + q) {1 - q^) {VW^ + l) I' ^ ^ ' 

Therefore, the symmetrized associator is characterized by one function Ts that is even in h, and has rational 
coefficients of h series expansion. We call this function Drinfeld prepotential, due to the role that the logarithm of 
the four point function plays in CFT, integrable systems and 4-dimensional SUSY models [T5] . 

6 Conclusion 

To conclude, we solved the regularized Knizhnik-Zamolodchikov equation and found the explicit expressions for 
associators in the case of the fundamental representation of gl{N). It turned out that associator is a finite series in 
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A and B with coefficients given by values of hypergeometric functions taken at the point z — 1. Also we presented 
some tests of our results. 

One possible direction of further research is to consider higher representations. For that case everything becomes 
more complicated due to more complicated relations for A and B. 

It would also be very interesting to understand the relation between conformal blocks and associator more 
deeply. 
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A Solution of KZ equation in the first two orders of perturbative 
expansion 

In this appendix we solve KZ equation in the first two orders of perturbative expansion. For KZ equation 

dG{z) _f^(A ^ B 



(^T^)=«. 



dz 

let us try to find the solution as a series in h: 

G{z) = G^{z) +hGi{z) + G2iz) + ... 
Equation ^ has a unique solution defined by boundary condition: 

Go(0) = 1 = 1(g) 1(g) 1 

We are interested in the value G(l). 
First order: 

1 1 

dGi{z) A B A f d.z ^ f dz 



0.il)=Af'^+Bf-^ 
dz z 1 — z J z J 1 — z 





The integrals here are divergent so we normalize them by: 

1 l-e 

-> 

e 

Then we have: 



Gi(z) = Alog(^) -i?log(^) 
and, keeping only finite and singular terms at point z = 1, 

Gi(l) = -Alog(e)-Blog(e) 

Second order: 

dG2(z) (A B\^,^ .9I, z ,1, 1-z ,l,z 1, 1-z 



. , , Gi(z) = = A^-log-- AB-log- +BA- log --B^- log- (57) 

dz \ z \ — z z e z 1 — e 1 — z e z 1 — e 
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Integrating over [e, z], we get: 

G2{z) = ]^A^ log' -^-AB (Li2(e) - Li2(z) ~ log(l - e) log + 

+ BA (Li2(e) - Li2(z) - log(l - z) log + ^B^ log' (58) 
At point z — 1 — e, for finite and singular terms we have: 

G2(l) = ^.42log'e + BAlogeloge+iB2log'e+[A,S]Li2(l) (59) 
Thus, up to the second order in h, G(l) has the following form (we drop the terms which tend to zero as e — >■ 0): 

G{z) = (^1 -;iSloge+ ^/I's'log'e^ {l + h^[A, B]U2{1)) (^1 - M log e + ^fi^A^ log' + O (h^) (60) 
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